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(57) For the purpose of enabling parallel imaging 
even when a navigator echo is used to phase-correct 
an imaging echo, the present invention involves: excit- 
ing spins within a subject to acquire an imaging echo 
generated by the excited spins along with a navigator 
echo, with a reduced field-of-view via a plurality of re- 
ceiver systems (802); conducting phase correction on 
the imaging echo based on the navigator echo (804); 
producing an intermediate image based on the phase- 
corrected imaging echo from each of the plurality of re- 
ceiver systems (806); generating a sensitivity matrix for 
the plurality of receiver systems (808); correcting the 
phase of matrix data in the sensitivity matrix (810); and 
producing an image with a full field-of-view based on the 
intermediate image and the phase-corrected sensitivity 
matrix (812). 
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Description 

[0001] The present invention relates to a magnetic 
resonance imaging apparatus, and more particularly to 
a magnetic resonance imaging apparatus that conducts 5 
parallel imaging. 

[0002] In a magnetic resonance imaging (MRI) appa- 
ratus, a subject to be imaged is carried into an internal 
space of a magnet system, i.e., an imaging space in 
which a static magnetic field is generated, a gradient 
magnetic field and a radio frequency (RF) magnetic field 
are applied to excite spins within the subject to thereby 
generate magnetic resonance signals, and an image is 
reconstructed based on the received signals. 
[0003] One scheme of the magnetic resonance imag- 
ing is parallel imaging. Parallel imaging is described by, 
for example, Klaas P. Pruessmann et al. in an article en- 
titled "SENSE: Sensitivity Encoding for Fast MRI" , Mag- 
netic Resonance in Medicine, 42: 952 — 962 (1999). 
[0004] Generally, in parallel imaging, magnetic reso- 
nance signals are acquired via a plurality of receiver sys- 
tems in a simultaneous and parallel manner. The acqui- 
sition of the magnetic resonance signals is conducted 
with a field-of-view (FOV) reduced by half, for example. 
By reducing the FOV by half, the rate of signal acquisi- 
tion is doubled. 

[0005] An image is reconstructed based on the sig- 
nals thus acquired. The image reconstruction is con- 
ducted in two steps. At the first step, an intermediate 
image is produced based on the signals acquired by the 
plurality of receiver systems. The image production em- 
ploys two-dimensional inverse Fourier transformation. 
The produced image has a reduced FOV. Because the 
FOV is reduced, aliasing images resulting from wrapa- 
round from outside of the FOV are contained in the im- 
age. 

[0006] At the second step, the aliasing images are 
brought back to their original positions by applying a cer- 
tain calculation to the image and an image with a whole 
FOV is produced. The calculation uses the following 
equation: 

[Equation 3] 
V = {S S)* 1 S* A, 

where 

V: pixel values of the image with the whole FOV, 
S: a sensitivity matrix, 
S*: an adjoint matrix of S, and 
A: pixel values of the intermediate image. 

[0007] The sensitivity matrix S is determined by the 
spatial distribution of the sensitivity of the plurality of re- 
ceiver systems. The sensitivity of a receiver system gen- 
erally has a complex form, and data of the sensitivity 



matrix therefore also has a complex form. Similarly, the 
pixel values V and A also have a complex form. 
[0008] One technique for taking a cross-sectional im- 
age of the head is MS-DW-EPI (multi-shot diffusion- 
weighted echo planar imaging). The technique involves 
capturing a cross-sectional image to which a weight is 
applied such that spins with less diffusion have a larger 
signal intensity, by a multi-shot (MS) echo planar imag- 
ing (EPI) technique. The captured cross-sectional im- 
age is suitable for diagnosing the presence of cerebral 
infarction. 

[0009] In the multi-shot echo planar imaging tech- 
nique, the collection of imaging echoes filling one screen 
image is spread over a plurality of times. Specifically, 
the spins are excited a plurality of number of times, and 
a certain number of imaging echoes are collected during 
each excitation. 

[0010] If positional shifting of the spins occurs be- 
tween the excitations due to pulsation of the brain in 
such a process, a phase error between the imaging ech- 
oes leads to ghosts generated in the reconstructed im- 
age. To avoid this, a navigator echo is collected for each 
excitation, and the phase of the navigator echo is used 
to correct the phase of the imaging echoes. 
[0011] It has been impossible to conduct parallel im- 
aging according to the MS-DW-EPI technique. This is 
because when the imaging echoes are phase-corrected 
by the navigator echo, the phase originating from the 
sensitivity of the receiver systems is unnecessarily cor- 
rected, resulting in inconsistency with the sensitivity ma- 
trix. 

[0012] The present invention therefore seeks to pro- 
vide a magnetic resonance imaging apparatus that en- 
ables parallel imaging even when a navigator echo is 
used to phase-correct an imaging echo. 
[0013] The present invention, for solving the afore- 
mentioned problem, is a magnetic resonance imaging 
apparatus characterized in comprising: acquiring 
means for exciting spins within a subject to acquire an 
imaging echo generated by the excited spins along with 
a navigator echo, with a reduced field-of-view via a plu- 
rality of receiver systems; first correcting means for con- 
ducting phase correction on said imaging echo based 
on said navigator echo; first image producing means for 
producing an intermediate image based on said phase- 
corrected imaging echo from each of said plurality of re- 
ceiver systems; generating means for generating a sen- 
sitivity matrix for said plurality of receiver systems; sec- 
ond correcting means for phase-correcting matrix data 
in said sensitivity matrix; and second image producing 
means for producing an image with a full field-of-view 
based on said intermediate image and said phase-cor- 
rected sensitivity matrix. 

[0014] Moreover, the invention for solving the afore- 
mentioned problem may be a magnetic resonance im- 
aging method characterized in comprising: exciting 
spins within a subject to acquire an imaging echo gen- 
erated by the excited spins along with a navigator echo, 
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with a reduced field-of-view via a plurality of receiver 
systems; conducting phase correction on said imaging 
echo based on said navigator echo; producing an inter- 
mediate image based on said phase-corrected imaging 
echo from each of said plurality of receiver systems; 5 
generating a sensitivity matrix for said plurality of receiv- 
er systems; phase-correcting matrix data in said sensi- 
tivity matrix; and producing an image with a full field-of- 
view based on said intermediate image and said phase- 
corrected sensitivity matrix. 

[0015] In the present invention, matrix data in a sen- 
sitivity matrix are phase-corrected, and therefore, incon- 
sistency with an imaging echo that is phase-corrected 
using a navigator echo is eliminated and an image can 
be correctly produced. 

[0016] The reduction factor for the reduced field-of- 
view preferably satisfies the following requirement so 
that a proper image may be obtained: 

[Equation 4] 
n £ R> 1 

where 

R: the reduction factor, and 
n: the number of receiver systems. 

[0017] The acquiring means preferably implements 
the reduced field-of-view by enlargement of sampling in- 
tervals for a fc-space so that the reduction of the field- 
of-view may be properly achieved. 
[0018] The acquiring means preferably implements 
the enlargement of the sampling intervals by enlarge- 
ment of a step difference of phase encoding when the 
field-of-view is to be reduced in the phase encoding di- 
rection. 

[0019] The plurality of receiver systems preferably 
have respective receiving coils so that the reception 
may be done in a simultaneous and parallel manner. 
[0020] The receiving coils are preferably surface coils 
so that the reception may be conducted close to the sub- 
ject. 

[0021] The acquiring means preferably employs an 
MS-DW-EPI technique in acquiring the imaging echo 
when an image incorporating diffusion is to be captured. 
[0022] The acquiring means preferably employs a 
technique other than the MS-DW-EPI technique in ac- 
quiring the imaging echo when an image other than that 
incorporating diffusion is to be captured. 
[0023] The generating means preferably generates 
the sensitivity matrix based on a spatial distribution of 
reception sensitivity of the plurality of receiver systems 
with respect to the full field-of-view so that the sensitivity 
matrix may be properly generated. 
[0024] The generating means preferably generates 
the sensitivity matrix after fitting the spatial distribution 
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of the magnitude of the reception sensitivity of each of 
the plurality of receiver systems to a two-dimensional 
polynomial so that the sensitivity matrix may be still 
more properly generated. 

[0025] The generating means preferably conducts the 
fitting by a method of least squares weighted depending 
upon the magnitude of the reception sensitivity so that 
the fitting may be properly achieved. 
[0026] A weight for the weighting is preferably the 
square of the magnitude of the reception sensitivity so 
that the fitting may be still more properly achieved. 
[0027] The two-dimensional polynomial is preferably 
a quadratic so that the fitting may be simplified. 
[0028] The second correcting means preferably ho- 
mogenizes the phase so that consistency with the phase 
correction on the imaging echo may be improved. The 
phase may be zero or a constant value other than zero. 
[0029] The second image producing means prefera- 
bly employs the following equation in producing the im- 
age so that an image with a full field-of-view may be 
properly obtained: 

[Equation 5] 

V = (S* S)" 1 S" A, 

where 

V: pixel values of the image with the full FOV, 

S: a sensitivity matrix, 

S*: an adjoint matrix of S, and 

A: pixel values of the intermediate image. 

[0030] Therefore, the present invention provides a 
magnetic resonance imaging apparatus that enables 
parallel imaging even when a navigator echo is used to 
phase-correct an imaging echo. 
[0031] Further objects and advantages of the present 
invention will be apparent from the following description 
of the preferred embodiments of the invention as illus- 
trated in the accompanying drawings. 
[0032] The invention will now be described in greater 
detail, by way of example, with reference to the draw- 
ings, in which:- 

Figure 1 is a block diagram of an apparatus in ac- 
cordance with one embodiment of the present in- 
vention. 

Figure 2 is a block diagram of an apparatus in ac- 
cordance with one embodiment of the present in- 
vention. 

Figure 3 shows a plurality of receiving coils. 
Figure 4 shows a plurality of receiver systems. 
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Figure 5 shows an exemplary pulse sequence exe- 
cuted by the apparatus in accordance with one em- 
bodiment of the present invention. 

Figure 6 shows a fc-space. 5 

Figure 7 shows a /c-space. 

Figure 8 is a flow chart of an operation of the appa- 
ratus in accordance with one embodiment of the 
present invention. 

Figure 9 is a flow chart of an operation of the appa- 
ratus in accordance with one embodiment of the 
present invention. 

Figure 10 is a functional block diagram of the appa- 
ratus in accordance with one embodiment of the 
present invention. 

Figure 11 is a functional block diagram of the appa- 
ratus in accordance with one embodiment of the 
present invention. 

Figure 12 is a functional block diagram of the appa- 
ratus in accordance with one embodiment of the 
present invention. 

Figure 13 shows an exemplary pulse sequence ex- 
ecuted by the apparatus in accordance with one 
embodiment of the present invention. 

[0033] Embodiments of the present invention will now 
be described in detail with reference to the accompany- 
ing drawings. Figure 1 shows a block diagram of a mag- 
netic resonance imaging apparatus, which is an embod- 
iment of the present invention. The configuration of the 
apparatus represents an embodiment of the apparatus 
in accordance with the present invention. 
[0034] As shown in Figure 1 , the present apparatus 
has a magnet system 1 00. The magnet system 1 00 has 
a main magnetic field coil section 102, a gradient coil 
section 106, and an RF coil section 108. These coil sec- 
tions have a generally cylindrical shape and are concen- 
trically disposed. A subject to be imaged 1 is rested on 
a cradle 500 and carried into and out of a generally cy- 
lindrical internal space (bore) of the magnet system 100 
by carrier means, which is not shown. 
[0035] The head of the subject 1 is fitted with a receiv- 
ing coil section 110. The receiving coil section 110 has 
a plurality of receiving coils. The receiving coil section 
110 will be described in detail later. 
[0036] The main magnetic field coil section 102 gen- 
erates a static magnetic field in the internal space of the 
magnet system 1 00. The direction of the static magnetic 
field is generally in parallel with the direction of the body 
axis of the subject 1 . That is, a magnetic field commonly 
referred to as a horizontal magnetic field is generated. 
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The main magnetic field coil section 102 is made using 
a superconductive coil, for example. However, the main 
magnetic field coil section 102 is not limited to the su- 
perconductive coil, but may be made using a normal 
conductive coil or the like. 

[0037] The gradient coil section 106 generates three 
gradient magnetic fields for imparting gradients to the 
static magnetic field strength in directions of three mu- 
tually perpendicular axes, i.e., a slice axis, a phase axis, 
and a frequency axis. 

[0038] When mutually perpendicular coordinate axes 
in the static magnetic field space are represented as X, 
Y, and Z, any one of the axes may be the slice axis. In 
this case, one of the two remaining axes is the phase 
axis and the other is the frequency axis. Moreover, the 
slice, phase, and frequency axes can be given arbitrary 
inclination with respect to the X-, Y-, and Z-axes while 
maintaining their mutual perpendicularity. In the present 
apparatus, the direction of the body axis of the subject 
1 is defined as the Z-axis direction. 
[0039] The gradient magnetic field in the slice axis di- 
rection is sometimes referred to as the slice gradient 
magnetic field. The gradient magnetic field in the phase 
axis direction is sometimes referred to as the phase en- 
coding gradient magnetic field. The gradient magnetic 
field in the frequency axis direction is sometimes re- 
ferred to as the readout gradient magnetic field. The re- 
adout gradient magnetic field is synonymous with the 
frequency encoding gradient magnetic field. In order to 
enable generation of such gradient magnetic fields, the 
gradient coil section 106 has three gradient coils, which 
are not shown. The gradient magnetic field will be some- 
times referred to simply as the gradient hereinbelow. 
[0040] The RF coil section 1 08 generates a radio fre- 
quency magnetic field in the static magnetic field space 
for exciting spins within the subject 1 . The generation of 
the radio frequency magnetic field will be sometimes re- 
ferred to as transmission of an RF excitation signal here- 
inbelow. Moreover, the RF excitation signal will be 
sometimes referred to as the RF pulse. 
[0041] Electromagnetic waves, i.e., magnetic reso- 
nance signals, generated by the excited spins are re- 
ceived by the receiving coil section 110. The magnetic 
resonance signals can also be received by the RF coil 
section 108. 

[0042] The magnetic resonance signals are those in 
a frequency domain, i.e., in a Fourier space. Since the 
magnetic resonance signals are encoded in two axes 
by the gradients in the phase- and frequency-axis direc- 
tions, the magnetic resonance signals are obtained as 
signals in a two-dimensional Fourier space. The phase 
encoding gradient and readout gradient are used to de- 
termine a position at which a signal is sampled in the 
two-dimensional Fourier space. The two-dimensional 
Fourier space will be sometimes referred to as the k- 
space hereinbelow. 

[0043] The gradient coil section 106 is connected with 
a gradient driving section 1 30. The gradient driving sec- 
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tion 1 30 supplies driving signals to the gradient coil sec- 
tion 106 to generate the gradient magnetic fields. The 
gradient driving section 130 has three driving circuits, 
which are not shown, corresponding to the three gradi- 
ent coils in the gradient coil section 106. 
[0044] The RF coil section 108 is connected with an 
RF driving section 140. The RF driving section 140 sup- 
plies driving signals to the RF coil section 1 08 to transmit 
the RF pulse, thereby exciting the spins within the sub- 
ject 1 . 

[0045] The receiving coil section 110 is connected 
with a data collecting section 1 50. The data collecting 
section 1 50 collects signals received by the receiving 
coil section 110 as digital data. The RF coil section 108 
can also be connected to the data collecting section 1 50 
to allow collection of signals received by the RF coil sec- 
tion 108. 

[0046] The gradient driving section 130, RF driving 
section 140 and data collecting section 150 are connect- 
ed with a sequence control section 160. The sequence 
control section 160 controls the gradient driving section 
130, RF driving section 140 and data collecting section 
1 50 to carry out the collection of magnetic resonance 
signals. 

[0047] The sequence control section 1 60 is, for exam- 
ple, constituted using a computer. The sequence control 
section 160 has a memory, which is not shown. The 
memory stores programs for the sequence control sec- 
tion 160 and several kinds of data. The function of the 
sequence control section 160 is implemented by the 
computer executing a program stored in the memory. 
[0048] The output of the data collecting section 1 50 is 
connected to a data processing section 170. Data col- 
lected by the data collecting section 1 50 are input to the 
data processing section 170. The data processing sec- 
tion 170 is, for example, constituted using a computer. 
The data processing section 170 has a memory, which 
is not shown. The memory stores programs for the data 
processing section 170 and several kinds of data. 
[0049] The data processing section 1 70 is connected 
to the sequence control section 160. The data process- 
ing section 170 is above the sequence control section 
160 and controls it. The function of the present appara- 
tus is implemented by the data processing section 170 
executing a program stored in the memory. 
[0050] The data processing section 1 70 stores the da- 
ta collected by the data collecting section 150 into the 
memory. A data space is established in the memory. The 
data space corresponds to the k-space. The data 
processing section 170 performs two-dimensional in- 
verse Fourier transformation on the data in the k-space 
to reconstruct an image. 

[0051] The data processing section 1 70 is connected 
with a display section 1 80 and an operating section 1 90. 
The display section 180 comprises a graphic display, 
etc. The operating section 190 comprises a keyboard, 
etc., provided with a pointing device. 
[0052] The display section 180 displays the recon- 



structed image output from the data processing section 
1 70 and several kinds of information. The operating sec- 
tion 1 90 is operated by a user, and the section 1 90 inputs 
several commands, information and so forth to the data 
5 processing section 1 70. The user interactively operates 
the present apparatus via the display section 180 and 
operating section 190. 

[0053] Figure 2 shows a block diagram of a magnetic 
resonance imaging apparatus of another type, which is 
one embodiment of the present invention. The configu- 
ration of the apparatus represents an embodiment of the 
apparatus in accordance with the present invention. 
[0054] The present apparatus has a magnet system 
1 00' of a type different from that of the apparatus shown 
in Figure 1 . Since the apparatus has a configuration sim- 
ilar to that of the apparatus shown in Figure 1 except for 
the magnet system 1 00', similar portions are designated 
by similar reference numerals and the explanation 
thereof will be omitted. 

[0055] The magnet system 1 00' has a main magnetic 
field magnet section 102\ a gradient coil section 106', 
and an RF coil section 108'. The main magnetic field 
magnet section 102' and the coil sections each consists 
of a pair of members facing each other across a space. 
These sections have a generally disk-like shape and are 
disposed to have a common center axis. A subject 1 is 
rested on a cradle 500 and carried into and out of an 
internal space (bore) of the magnet system 100' by car- 
rier means, which is not shown. 
[0056] The head of the subject 1 is fitted with a receiv- 
ing coil section 110. The receiving coil section 110 has 
a plurality of receiving coils. The receiving coil section 
110 will be described in detail later. 
[0057] The main magnetic field magnet section 102' 
generates a static magnetic field in the internal space of 
the magnet system 1 00*. The direction of the static mag- 
netic field is generally orthogonal to the direction of the 
body axis of the subject 1 . That is, a magnetic field com- 
monly referred to as a vertical magnetic field is gener- 
ated. The main magnetic field magnet section 102' is 
made using a permanent magnet, for example. Howev- 
er, the main magnetic field magnet section 102' is not 
limited to a permanent magnet, but may be made using 
a super or normal conductive electromagnet or the like. 
[0058] The gradient coil section 106' generates three 
gradient magnetic fields for imparting gradients to the 
static magnetic field strength in directions of three mu- 
tually perpendicular axes, i.e., a slice axis, a phase axis 
and a frequency axis. 

[0059] When mutually perpendicular coordinate axes 
in the static magnetic field space are represented as X, 
Y, and Z, any one of the axes may be the slice axis. In 
this case, one of the two remaining axes is the phase 
axis and the other is the frequency axis. Moreover, the 
slice, phase, and frequency axes can be given arbitrary 
inclination with respect to the X-, Y-, and Z-axes while 
maintaining their mutual perpendicularity. In the present 
apparatus, the direction of the body axis of the subject 
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1 is again defined as the Z-axis direction. 
[0060] The gradient magnetic field in the slice axis di- 
rection is sometimes referred to as the slice gradient 
magnetic field. The gradient magnetic field in the phase 
axis direction is sometimes referred to as the phase en- 5 
coding gradient magnetic field. The gradient magnetic 
field in the frequency axis direction is sometimes re- 
ferred to as the readout gradient magnetic field. The re- 
adout gradient magnetic field is synonymous with the 
frequency encoding gradient magnetic field. In order to 
enable generation of such gradient magnetic fields, the 
gradient coil section 106' has three gradient coils, which 
are not shown. 

[0061] The RF coil section 1 08' transmits an RF pulse 
to the static magnetic field space for exciting spins within 
the subject 1. Electromagnetic waves, i.e., magnetic 
resonance signals, generated by the excited spins are 
received by the receiving coil section 110. The magnetic 
resonance signals can also be received by the RF coil 
section 108'. The signals received by the receiving coil 
section 110 or RF coil section 108' are input to the data 
collecting section 150. 

[0062] The receiving coil section 110 will now be de- 
scribed. The receiving coil section 110 has a plurality of 
receiving coils 112 — 118 disposed close to the head of 
the subject 1 as shown in Figure 3. Each receiving coil 
112 — 118 forms a closed loop. The receiving coils 
112 — 118 are disposed surrounding the head of the 
subject 1 . For the receiving coils 1 1 2 — 1 1 8, surface coils 
are employed, for example. The surface coils are suita- 
ble for receiving signals very close to the subject. 
[0063] The receiving coils 112 — 118 are configured 
so that mutual interference is substantially prevented. A 
plurality of receiving coils having no mutual interference 
are sometimes referred to collectively as a phased-array 
coil. The number of receiving coils constituting the 
phased-array coil may be any plural number greater or 
smaller than four. The following description will be made 
on a case in which the number of receiving coils is four 
but also applies to a case in which the number of receiv- 
ing coils is a plural number greater or smaller than four. 
[0064] The receiving coils 112 — 118 constituting a 
phased-array coil separately receive respective mag- 
netic resonance signals. The signals received by the re- 
ceiving coils 112 — 118 are input to receiver circuits 
152 — 158 in the data collecting section 150 as shown 
in Figure 4. 

[0065] The receiving coil 1 1 2 and receiver circuit 1 52 
constitute a receiver system 1. The receiving coil 114 
and receiver circuit 154 constitute a receiver system 2. 
The receiving coil 11 6 and receiver circuit 156 constitute 
a receiver system 3. The receiving coil 118 and receiver 
circuit 158 constitute a receiver system 4. 
[0066] By such a plurality of receiver systems having 
separate receiving coils, signals can be received by the 
plurality of receiver systems in a simultaneous and par- 
allel manner. The receiver systems 1 — 4 represent an 
embodiment of the plurality of receiver systems in ac- 



cordance with the present invention. 
[0067] Figure 5 shows a pulse sequence of a scan ac- 
cording to the MS-DW-EPI technique. The pulse se- 
quence proceeds from the left to the right. The same 
applies in the following description. In Figure 5, (1) 
shows a pulse sequence of RF signals. (2) — (5) show 
pulse sequences of gradient magnetic fields. (2) repre- 
sents a slice gradient, (3) represents a frequency en- 
coding gradient, (4) represents a phase encoding gra- 
dient, and (5) represents a motion probing gradient. It 
should be noted that the static magnetic field is con- 
stantly applied at a fixed magnetic field strength. 
[0068] First, spin excitation is effected by a 90° pulse. 
After a certain time period from the 90° excitation, 180° 
excitation is effected by a 180° pulse. These are selec- 
tive excitations under a slice gradient G sllce . Before and 
after the 1 80° excitation, a motion probing gradient G mpg 
is applied. This applies weighting depending upon dif- 
fusion of spins to a magnetic resonance signal. 
[0069] Next, a frequency encoding gradient G freq is 
solely applied to read a navigator echo. A pair of the 
navigator echoes are read. Each navigator echo is rep- 
resented by its central echo. The same applies in the 
following description. 

[0070] Next, a phase encoding gradient G phase and a 
frequency encoding gradient G freq are applied in a pre- 
defined sequence and a plurality of imaging echoes are 
sequentially read. The plurality of imaging echoes have 
different phase encodings. Again, each imaging echo is 
represented by its central echo. The same applies in the 
following description. 

[0071] Such a pulse sequence is repeated a prede- 
termined number of times in a cycle of a repetition time 
TR, and a pair of navigator echoes and a plurality of im- 
aging echoes are read for each repetition. That is, a mul- 
ti-shot scan is performed. The phase encoding for the 
imaging echo is changed for each repetition, and imag- 
ing echoes filling one screen image are acquired in the 
predetermined number of times of repetitions. 
[0072] The imaging echoes in each repetition are 
phase-collected by a navigator echo in the same repe- 
tition. This removes any phase error due to pulsation of 
the brain from the imaging echoes. 
[0073] By reading the imaging echoes with the phase 
and frequency encodings, data in a /c-space are sam- 
pled. Figure 6 shows a conceptual diagram of the /c- 
space. As shown, the horizontal axis kx of the /c-space 
is a frequency axis and the vertical axis ky is a phase 
axis. 

[0074] In Figure 6, each of a plurality of laterally ex- 
tending rectangles represents a data sampling position 
on the phase axis. The number shown in each rectangle 
represents the amount of phase encoding. The amounts 
of phase encoding are normalized by tUN. N is the 
number of times of sampling in the phase axis direction. 
[0075] The amount of phase encoding is zero at the 
center of the phase axis ky. The amount of phase en- 
coding increases from the center to both ends. The po- 
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larities of the increases are opposite to each other. The 
sampling interval, i.e., the step difference between the 
amounts of phase encoding is n/N. By performing two- 
dimensional inverse Fourier transformation on the data 
in such a /c-space, a cross-sectional image is recon- 5 
structed. The reconstructed image represents a whole 
field-of-view (FOV). The whole FOV will be sometimes 
referred to as the full FOV hereinbelow. 
[0076] In parallel imaging, the sampling intervals for 
the /c-space are increased to reduce the number of times 
of sampling for improving the speed of imaging. Specif- 
ically, sampling in the phase axis ky direction is conduct- 
ed in an 

[0077] alternate manner, for example, to reduce the 
number of times of sampling by half as indicated by 
hatching shown in Figure 7. This reduces the imaging 
time by half, thereby improving the speed of imaging. 
[0078] By the alternate sampling, the sampling inter- 
vals are doubled. By doubling the sampling intervals, the 
FOV of the reconstructed image is reduced by half as 
compared with the full FOV 

[0079] The doubling of the sampling intervals in the 
phase encoding direction is achieved by setting the step 
difference between the amounts of phase encoding to 
2idN. This reduces the FOV by half in the phase encod- 
ing direction. 

[0080] Generally, when the sampling interval, or the 
step difference between the amounts of phase encod- 
ing, is increased by a factor of R, the FOV is reduced by 
a factor of R. R is sometimes referred to as a reduction 
factor. In Figure 7, R = 2. 

[0081 ] When the number of receiver systems is n, the 
reduction factor R preferably satisfies the following re- 
lationship so that an output image with a full FOV as 
described later may be properly obtained: 

[Equation 6] 
n Z R > 1, 

where 

R: the reduction factor, and 
n: the number of receiver systems. 

[0082] The operation of the present apparatus will 
now be described. Figure 8 shows a flow chart of the 
operation of the present apparatus. As shown, at Step 
701 , reception sensitivity distribution measurement is 
conducted. At this step, the spatial distribution of the 
sensitivity of the plurality of receiver systems is meas- 
ured. 

[0083] The spatial distribution of the sensitivity of the 
receiver systems is obtained as sensitivity map images. 
The sensitivity map images are produced from images 
obtained by scanning the same slice through the subject 
1 using the RF coil section 108 (108') and receiving coil 



section 110, for example. 

[0084] Specifically, an image captured using the RF 
coil section 108 (108') is defined as a reference, images 
captured using the receiving coils 112 — 118 are de- 
fined as measured images, and the sensitivity maps are 
generated by, for example, calculating the ratio of the 
measured images and reference image on a pixel-by- 
pixel basis. The reference image and measured images 
are captured by scanning of a full FOV. Thus, the sen- 
sitivity map images are obtained for the receiving coils 
112 — 118 with a full FOV. Such scanning is sometimes 
referred to as a calibration scan. 
[0085] Next, at Step 703, sensitivity matrix generation 
is conducted. The sensitivity matrix is generated based 
on the sensitivity map image of each receiving coil. 
Since the sensitivity map images are obtained with the 
full FOV, the sensitivity matrix can be properly generat- 
ed. The sensitivity map image will be sometimes re- 
ferred to simply as the sensitivity map hereinbelow. 
[0086] The sensitivity matrix is an n x R matrix. The 
symbol n designates the number of receiver systems, 
and R designates the reduction factor. When n = 4 and 
R = 2, the sensitivity matrix S is as follows: 



[Equation 7] 



S = 



[0087] In the sensitivity matrix S, s11, s21, s31 and 
s41 designate values of the same pixel in the sensitivity 
map images of the receiving coils 112,114,116 and 1 1 8. 
Values of a pixel in the sensitivity maps that lies at a 
distance of 1/2 FOV from the former pixel in the phase 
encoding direction are designated by s12, s22, s32 and 
s42. These values are complex. 
[0088] Next, at Step 705, sensitivity matrix phase cor- 
rection is conducted. The phase correction is achieved 
by setting the phase of complex data s in the sensitivity 
matrix S to zero or a predetermined constant value. Spe- 
cifically, when the complex data s is defined as: 

[Equation 8] 
s = Me^, 
the phase is defined as follows: 

[Equation 9] 
9=0, 
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or 

[Equation 10] 
6 = const (* 0). 

[0089] Next, at Step 707, a scan is conducted. The 
scan is conducted according to the MS-DW-EPI tech- 
nique. The scan according to the MS-DW-EPI technique 
is conducted on a reduced FOV by increasing the sam- 
pling intervals for the fr-space. The reduced FOV is, for 
example, a 1/2 FOV. The reduction factor R is not limited 
to 1/2 but may be any appropriate value. The navigator 
echo and imaging echoes are received via the plurality 
of receiver systems 1 — 4 in a simultaneous and parallel 
manner. 

[0090] Next, at Step 709, phase correction on the im- 
aging echoes is conducted. The phase correction is 
achieved by phase-correcting the imaging echoes by 
the navigator echo for each shot of multi shots. 
[0091] Next, at Step 711 , intermediate image produc- 
tion is conducted. An intermediate image is produced 
by performing two-dimensional inverse Fourier transfor- 
mation on the phase-corrected imaging echoes of the 
plurality of receiver systems. The intermediate image 
contains aliasing images because it has a reduced FOV. 
[0092] Next, at Step 713, output image production is 
conducted. The output image is produced by a calcula- 
tion using the intermediate image and sensitivity matrix. 
The output image production uses the following equa- 
tion. The equation is similar to that described in the ar- 
ticle cited earlier. 

[Equation 11] 

V = (S* S)" 1 S* A, 

where 

V: pixel values of the image with a full FOV, 

S: a sensitivity matrix, 

S*: an adjoint matrix of S, and 

A: pixel values of the intermediate image. 

[0093] Although the pixel values A in the intermediate 
image are values after phase-correction in this equation, 
the calculation can be carried out without inconsistency 
because the sensitivity matrix S is also phase-corrected. 
Thus, a cross-sectional image in which the aliasing im- 
ages are rearranged to their original positions can be 
obtained with a full FOV. 

[0094] The cross-sectional image is displayed and 
stored at Step 715. The display of the cross-sectional 
image is executed by the display section 180, and the 
storage thereof is directed to the memory in the data 
processing section 170. 



[0095] Moreover, the phase correction may be con- 
ducted on the sensitivity maps because the data in the 
sensitivity matrix are exactly the data of the sensitivity 
maps. The flow chart in this case is shown in Figure 9. 
5 Figure 9 is different from that shown in Figure 8 only in 
Step 705*. 

[0096] As shown in Figure 9, the phase correction is 
conducted on the sensitivity maps at Step 705', and the 
sensitivity matrix is generated based on the phase-cor- 
10 rected sensitivity maps at Step 703. The sensitivity ma- 
trix thus generated has the corrected phase. That is, the 
phase correction on the sensitivity matrix is achieved in 
effect. 

[0097] Figure 1 0 shows a functional block diagram of 

15 the present apparatus for the above-described opera- 
tion. As shown, the apparatus comprises a signal ac- 
quiring section 802, a phase correcting section 804, an 
intermediate image producing section 806, a sensitivity 
matrix generating section 808, a phase correcting sec- 

20 tion 810 and an output image producing section 812. 
[0098] The signal acquiring section 802 acquires nav- 
igator echoes and imaging echoes. The echoes are ac- 
quired according to the MS-DW-EPI technique with a 
reduced FOV. The signal acquiring section 802 corre- 

25 sponds to the function of the present apparatus de- 
scribed at Step 707 shown in Figure 8. The signal ac- 
quiring section 802 is an embodiment of the acquiring 
means of the present invention. 
[0099] The phase correcting section 804 conducts the 

30 phase correction on the imaging echoes. The phase cor- 
recting section 804 corresponds to the function of the 
present apparatus described at Step 709 shown in Fig- 
ure 8. The phase correcting section 804 is an embodi- 
ment of the first correcting means of the present inven- 

35 tion. 

[0100] The intermediate image producing section 806 
conducts the intermediate image production based on 
the phase-corrected imaging echoes. The intermediate 
image producing section 806 corresponds to the func- 

40 tion of the present apparatus described at Step 711 
shown in Figure 8. The intermediate image producing 
section 806 is an embodiment of the first image produc- 
ing means of the present invention. 
[0101] The sensitivity matrix generating section 808 

45 conducts the sensitivity matrix generation. The sensitiv- 
ity matrix generating section 808 corresponds to the 
functions of the present apparatus described at Steps 
701 and 703 shown in Figure 8. The sensitivity matrix 
generating section 808 is an embodiment of the gener- 

50 ating means of the present invention. 

[0102] The sensitivity matrix generating section 808 
is comprised of a reception sensitivity distribution meas- 
uring section 902 and a matrix generating section 904 
as shown in Figure 11. The reception sensitivity distri- 

55 bution measuring section 902 corresponds to the func- 
tion of the present apparatus described at Step 701 
shown in Figure 8. The matrix generating section 904 
corresponds to the function of the present apparatus de- 
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scribed at Step 703 shown in Figure 8. 
[0103] The phase correcting section 810 conducts the 
phase correction on the sensitivity matrix. The phase 
correcting section 810 corresponds to the function of the 
present apparatus described at Step 705 shown in Fig- 
ure 8 or Step 705' shown in Figure 9. The phase cor- 
recting section 81 0 is an embodiment of the second cor- 
recting means of the present invention. 
[0104] The output image producing section 81 2 con- 
ducts the output image production using the intermedi- 
ate image and sensitivity matrix. The output image pro- 
ducing section 812 corresponds to the function of the 
present apparatus described at Step 713 shown in Fig- 
ure 8. The output image producing section 812 is an em- 
bodiment of the second image producing means of the 
present invention. 

[0105] The sensitivity matrix generating section 808 
may comprise a fitting section 906 between the recep- 
tion sensitivity distribution measuring section 902 and 
matrix generating section 904 as shown in Figure 12. 
[0106] The fitting section 906 conducts functional fit- 
ting on the magnitude in the sensitivity maps obtained 
by the reception sensitivity distribution measuring sec- 
tion 902. In the fitting, the phases are preferably correct- 
ed beforehand so that the fitting may be facilitated. This 
also achieves the phase correction on the sensitivity 
matrix. 

[01 07] The function used in the fitting is a two-dimen- 
sional polynomial. The polynomial is preferably quad- 
ratic so that the order is small and the fitting is simplified. 
However, the order may be high, such as an order of 
three or greater. 

[0108] The fitting in the case of an order of two is 
achieved by calculating coefficients a — f of the follow- 
ing two-dimensional quadratic polynomial by a method 
of least squares: 

[Equation 12] 
M = ax 2 + bx 2 + cx + dx + exy + f 

[0109] At that time, the square error is preferably 
weighted according to the magnitude of the sensitivity 
maps so that an effect of noise on the fitting is mitigated. 
Fitting with higher accuracy can be achieved by setting 
the weight to the square of the magnitude. 
[0110] Since the sensitivity maps can be mathemati- 
cally expressed by the fitting, the output image can be 
produced using the sensitivity matrix S without a large 
error even if positional shifting of the pixels occurs be- 
tween the calibration scan image (sensitivity maps) and 
the actually scanned image (intermediate image) due to 
the motion of the subject 1 or the like. That is, parallel 
imaging less affected by motion can be achieved. 
[01 1 1] The reason of this is that, by using the mathe- 
matical expression, it is possible to continuously extrap- 
olate the sensitivity value not only for an exact pixel po- 



sition but for proximate pixel positions, and therefore, a 
sensitivity value having high correctness can be ob- 
tained even if the pixel position is shifted. 
[0112] Moreover, in general, fitting including the 

5 phase is difficult because receiver systems (coils) in a 
phased-array coil have different phase characteristics; 
however, the phase of the sensitivity matrix is homoge- 
nized in the present invention and the phase may be 
regarded as having been fitted in a sense, and therefore, 

10 exact fitting of the sensitivity matrix is possible. 

[0113] Parallel imaging involving navigator echo col- 
lection is not limited to the MS-DW-EPI. Parallel imaging 
involving navigator echo collection may include, besides 
the MS-DW-EPI, imaging according to a three-dimen- 

15 sional gradient echo technique, for example. 

[0114] Figure 1 3 shows a pulse sequence of the tech- 
nique. In Figure 13, (1) shows a sequence of an RF 
pulse. (2) — (4) show pulse sequences of gradient mag- 
netic fields. (2) represents a slice gradient and a phase 

20 encoding gradient in the slice direction, (3) represents 
a frequency encoding gradient, and (4) represents a 
phase encoding gradient. It should be noted that the 
static magnetic field is constantly applied at a fixed mag- 
netic field strength. 

25 [0115] First, spin excitation is effected by an a 0 pulse. 
The a° excitation is selective excitation under a slice 
gradient G S | ice . After the a° excitation, a phase encoding 
gradient G S | ice in the slice direction, a frequency encod- 
ing gradient G freq and a phase encoding gradient G phase 

30 are applied in a predefined sequence and an imaging 
echo is read. 

[0116] After reading the imaging echo, the phase en- 
coding gradient G s[jce in the slice direction and phase 
encoding gradient G phase are rewound. Thereafter, a 
35 frequency encoding gradient G freq is solely applied to 
read a navigator echo. 

[0117] Such a pulse sequence is repeated a prede- 
termined number of times in a cycle of a repetition time 
TR, and an imaging echo and a navigator echo are read 

40 for each repetition. The phase encoding for the imaging 
echo is changed for each repetition, and imaging ech- 
oes filling one screen image are acquired in the prede- 
termined number of repetitions. The imaging echo in 
each repetition is phase-collected by a navigator echo 

45 jn the same repetition. 

[0118] By reading the imaging echo using the phase 
encodings in two directions and the frequency encoding, 
data in a three-dimensional /c-space are sampled. By in- 
creasing the sampling intervals, the data are collected 

50 with a three-dimensional reduced FOV. 

[01 1 9] By performing three-dimensional inverse Fou- 
rier transformation on the data, a three-dimensional im- 
age is reconstructed. The three-dimensional image is 
an intermediate image with the reduced FOV. From the 

55 intermediate image, an output image is produced with 
a full FOV using a sensitivity matrix S. It should be noted 
that the sensitivity matrix used is of a three-dimensional 
type. 
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[0120] In imaging the abdomen, body motion be- 
comes an issue that makes a calibration scan and an 
actual scan susceptible to positional shifting; however, 
according to the present invention, parallel imaging can 
be achieved unaffected by such positional shifting. 
[0121] For the sake of good order, various aspects of 
the invention are set out in the following clauses:- 

1 . A magnetic resonance imaging apparatus com- 
prising: 

an acquiring device for exciting spins within a 
subject to acquire an imaging echo generated 
by the excited spins along with a navigator 
echo, with a reduced field-of-view via a plurality 
of receiver systems; 

the first correcting device for conducting phase 
correction on said imaging echo based on said 
navigator echo; 

the first image producing device for producing 
an intermediate image based on said phase- 
corrected imaging echo from each of said plu- 
rality of receiver systems; 
a generating device for generating a sensitivity 
matrix for said plurality of receiver systems; 
the second correcting device for phase-correct- 
ing matrix data in said sensitivity matrix; and 
the second image producing device for produc- 
ing an image with a full field-of-view based on 
said intermediate image and said phase-cor- 
rected sensitivity matrix. 

2. The magnetic resonance imaging apparatus of 
clause 1 , wherein 

the reduction factor for said reduced field-of- 
view satisfies the following requirement: 

n £ R> 1, 

where 

R: the reduction factor, and 

n: the number of receiver systems. 

3. The magnetic resonance imaging apparatus of 
clause 1 , wherein said acquiring device implements 
said reduced field-of-view by enlargement of sam- 
pling intervals for a /c-space. 

4. The magnetic resonance imaging apparatus of 
clause 3, wherein said acquiring device implements 
said enlargement of the sampling intervals by en- 
largement of a step difference of phase encoding. 

5. The magnetic resonance imaging apparatus of 
clause 1, wherein said plurality of receiver systems 
have respective receiving coils. 



6. The magnetic resonance imaging apparatus of 
clause 5, wherein said receiving coils are surface 
coils. 

5 7. The magnetic resonance imaging apparatus of 
clause 1 , wherein said acquiring device employs an 
MS-DW-EPI technique in acquiring said imaging 
echo. 

10 8. The magnetic resonance imaging apparatus of 
clause 1, wherein said acquiring device employs a 
technique other than the MS-DW-EPI technique in 
acquiring said imaging echo. 

15 9. The magnetic resonance imaging apparatus of 
clause 1 , wherein said generating device generates 
said sensitivity matrix based on a spatial distribution 
of reception sensitivity of said plurality of receiver 
systems with respect to the full field-of-view. 

20 

10. The magnetic resonance imaging apparatus of 
clause 9, wherein said generating device generates 
said sensitivity matrix after fitting the spatial distri- 
bution of the magnitude of the reception sensitivity 

25 of each of said plurality of receiver systems to a two- 
dimensional polynomial. 

11. The magnetic resonance imaging apparatus of 
clause 10, wherein said generating device conducts 

30 said fitting by a method of least squares weighted 
depending upon the magnitude of the reception 
sensitivity. 

12. The magnetic resonance imaging apparatus of 
35 clause 1 1 , wherein a weight for said weighting is the 

square of the magnitude of the reception sensitivity. 

13. The magnetic resonance imaging apparatus of 
clause 10, wherein said two-dimensional polynomi- 

<o al is a quadratic. 

14. The magnetic resonance imaging apparatus of 
clause 1 , wherein said second correcting device ho- 
mogenizes said phase. 

45 

15. The magnetic resonance imaging apparatus of 
clause 1, wherein said second correcting device 
sets said phase to zero. 

50 16. The magnetic resonance imaging apparatus of 
clause 1, wherein said second correcting device 
sets said phase to a constant value other than zero. 

17. The magnetic resonance imaging apparatus of 
55 clause 1 , wherein said second image producing de- 
vice employs the following equation in producing 
said image: 
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V = (S' S)* 1 S* A, 

where 

V: pixel values of the image with the full FOV, 

S: a sensitivity matrix, 

S*: an adjoint matrix of S, and 

A: pixel values of the intermediate image. 



Claims 

1 . A magnetic resonance imaging apparatus compris- 
ing: 

an acquiring device for exciting spins within a 
subject to acquire an imaging echo generated 
by the excited spins along with a navigator 
echo, with a reduced field-of-view via a plurality 
of receiver systems; 

the first correcting device for conducting phase 
correction on said imaging echo based on said 
navigator echo; 

the first image producing device for producing 
an intermediate image based on said phase- 
corrected imaging echo from each of said plu- 
rality of receiver systems; 
a generating device for generating a sensitivity 
matrix for said plurality of receiver systems; 
the second correcting device for phase-correct- 
ing matrix data in said sensitivity matrix; and 
the second image producing device for produc- 
ing an image with a full field-of-view based on 
said intermediate image and said phase-cor- 
rected sensitivity matrix. 



largement of a step difference of phase encoding. 

5. The magnetic resonance imaging apparatus of any 
preceding claim, wherein said plurality of receiver 

5 systems have respective receiving coils. 

6. The magnetic resonance imaging apparatus of 
claim 5, wherein said receiving coils are surface 
coils. 

10 

7. The magnetic resonance imaging apparatus of any 
preceding claim, wherein said acquiring device em- 
ploys an MS-DW-EPI technique in acquiring said 
imaging echo. 

15 

8. The magnetic resonance imaging apparatus of any 
one of claims 1 to 6, wherein said acquiring device 
employs a technique other than the MS-DW-EPI 
technique in acquiring said imaging echo. 

20 

9. The magnetic resonance imaging apparatus of any 
preceding claim, wherein said generating device 
generates said sensitivity matrix based on a spatial 
distribution of reception sensitivity of said plurality 

25 of receiver systems with respect to the full field-of- 
view. 

10. The magnetic resonance imaging apparatus of 
claim 9, wherein said generating device generates 

30 said sensitivity matrix after fitting the spatial distri- 
bution of the magnitude of the reception sensitivity 
of each of said plurality of receiver systems to a two- 
dimensional polynomial. 

35 



2. The magnetic resonance imaging apparatus of 
claim 1 , wherein 

the reduction factor for said reduced field-of- 40 
view satisfies the following requirement: 



n £ R > 1, 

where 



45 



R: the reduction factor, and 

n: the number of receiver systems. 

3. The magnetic resonance imaging apparatus of 
claim 1 or 2, wherein said acquiring device imple- 
ments said reduced field-of-view by enlargement of 
sampling intervals for a /c-space. 

4. The magnetic resonance imaging apparatus of 
claim 3, wherein said acquiring device implements 
said enlargement of the sampling intervals by en- 
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